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ABSTRACT 

After insertion into the near-earth parking orbit a Go, No-Go decision has 
to be made, based on one minute of ship's tracking data. The deciding factor 
for the Go, No-Go decision is the perigee height of the parking orbit; a mini- 
mum perigee height of 75 n.mi. is required. In this paper, a statistical analysis 
has been made of the error in the perigee height, as determined from one min- 
ute of ship's tracking data. It is shown, that the computed perigee height has to 
be 0.5 to 17 n.mi. higher than 75 n.mi. in order to insure that the actual perigee 
height exceeds 75 n.mi. with 99.5% probability. The use of a 17.n.mi. padding 
is recommended in order to cover worst case conditions. 
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SUMMARY 


At insertion, the following parameters: 

radius r 

O 

speed v o 

flight path angle y o 

will be determined from one minute of a ship's tracking data. In reference 2, it 
is required that these parameters be determined to the following accuracies: 


3 CJ r 0 

- 2.4 n.mi . 

(4.44 km) 

3 CT v 0 

= 16 ft/sec 

(4. 87 m/sec) 

3 cr 

7 O 

= 0.16° 

(2.79 mrad) 


A statistical analysis of the perigee height error due to the above insertion 
parameter uncertainties is presented in this paper. Gaussian distributions are 
assumed for the errors in the insertion parameters (i.e., measured minus actual) 
and correlations between the errors are taken into account. The results are ap- 
plied to the Apollo Go, No-Go decision at insertion. 

For a Go decision, an actual perigee height of 75 n.mi. is required with 
99.5% probability (reference 1). It is shown that the computed perigee height 
(from ship's tracking data) has to be 0.5 to 17 n.mi. higher than 75 n.mi. in 
order to meet this requirement. The required padding is a function of the ec- 
centricity e of the parking orbit, the true anomaly 6 at insertion, the errors in 
the insertion parameters and their correlations. 

The results of the analysis may be summed up in the following recommen- 
dations: 


1. A padding of 17 n.mi. should be used to insure with a 99.5% probability 
that a Go decision is correct. This value covers the following cases: 

a. All eccentricities in the region 0 < e £ 0.01. 

b. All values of true anomaly at insertion. 

c. All values of correlation coefficients p between errors of the in- 
sertion parameters in the region 0 < | p | < 0.9. 

2. The parking orbit should have the smallest possible eccentricity and a 
true anomaly close to 0° at insertion in order to minimize the perigee 
height error. 
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THE EFFECT OF INSERTION SHIP 
TRACKING ERRORS ON THE APOLLO 
GO, NO-GO DECISION 


1. INTRODUCTION 

For the Apollo mission a decision to continue or not continue, i.e., the Go, 
No-Go decision, will have to be made shortly after the spacecraft is inserted 
into the parking orbit. 

Three sources of data are available for the Go, No-Go decision: the inser- 
tion ship and two on board inertial guidance systems. One of the on board sys- 
tems is located in the S-IVB and the other in the Apollo spacecraft (Command 
Module). Normally, the Go, No-Go decision will be made using the two on board 
systems. In the event these two systems are in disagreement, the insertion 
ship will be used as an arbiter (reference 1). 

This report is a statistical analysis of the effect of tracking measurement 
errors by the insertion ship on the error in perigee for near-earth Apollo 
parking orbits. The results of this analysis are applied to the Apollo Go, No-Go 
decision at insertion. 

Figure 1 shows the geometry at insertion. The critical insertion parameters 
are the radius r o , the magnitude of the velocity v o , and the flight path angle y . 

As required in reference 2, these insertion parameters must be determined 
from the ship's tracking data with the following accuracies: 



- 2.4 n.mi 

(4. 44 km) 

3CT v 0 

= 16 ft/sec 

(4.87 m/ sec) 

3 "r 0 

= 0.16° 

(2. 79 mrad) 


It is shown in reference 3, Chapter 3, that these accuracy requirements 
can be met with one minute of tracking by the ship if the tracking error model 
from reference 4 is used. Three minutes of tracking is obtainable with the 
positions of the insertion ship as planned in reference 3. Thus, since only one 
minute of tracking data is needed, this provides a margin for loss of data due 
to delayed acquisition, interrupted tracking, etc. 

The Go, No-Go decision will be based upon the perigee height. In reference 
2 a minimum perigee height of 71 n.mi. (131 km) is required. In order to allow 
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INSERTION 



Figure 1. Insertion Geometry. 

for two orbits before re-entry, a minimum perigee height of 75 n.mi. (139 km) 
is suggested in reference 1, and a 99.5% probability is required to achieve this 
minimum height if a Go decision is made (reference 1). 

The perigee height may be calculated from the insertion parameters by 
means of the Keplerian two-body equations of motion. However, due to errors 
in the insertion parameters as determined from the ship's tracking data, there 
will be an error associated with the calculated perigee height. It is the purpose 
of this report to study the probability distribution of perigee error and apply the 
results to the Apollo Go, No-Go decision at insertion. 

Distribution determination of the perigee error is important from the stand- 
point of astronaut safety. If the actual perigee height is lower than the calculated 
perigee height , an incorrect Go decision could be made (i.e., one where the 
actual perigee height is lower than 75 n.mi.). This would be detrimental to the 
safety of the astronauts due to the excessive re-entry heat. 
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As stated in reference 3, a circular parking orbit is planned. However, an 
elliptical orbit should not be ruled out, in which case the perigee height may be 
at 85 n.mi. (157 km) and the apogee height at 150 n.mi. (278 km). Both circular 
and elliptical parking orbits are therefore treated in this analysis. 

2. THE CUMULATIVE DISTRIBUTION OF PERIGEE ERROR 

2.1 The Perigee Error 


The spacecraft has the actual or true insertion parameters r o , v o , and y . 
The actual perigee radius r p can be determined from the Keplerian two -body 
equations of motion. 
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2 

O 



( 2 . 1 . 1 ) 


(1 - e 2 ) 



( 2 . 1 . 2 ) 


r p - a (1 - e) 


(2.1.3) 


where a = semi-major axis of the parking orbit 

e = eccentricity of the parking orbit 

ix = gravitational constant of the earth (398603.2 (km) 3 / sec 2 ) 


by eliminating a and e as shown in Appendix A. Thus, the perigee radius r 
is a function of r , v , and y . 

o o ' o 


Since the insertion parameters are determined from the ship's tracking 
data, they will be in error by the amounts Ar o , Av o , and Ay o . Therefore, the 
calculated perigee radius r p ^ will deviate by Ar p from the actual perigee 
radius r . 

p 


- r 


Pc a 1 


r (r + Ar , v + Av , y + Ay ) - r (r , v , y ) 

p v o o’ o o / o / o / p v o’ o’ / o / 


(2.1.4) 

This equation determines A r as a function of the insertion parameters and 
their errors. By eliminating r , v , and y with the aid of equations (2.1.1) 
through (2.1.3), an alternate expression for Ar results (equations B.63, B.64, 
and B.68 of Appendix B). 

Ar p - A r p (a, e, 6 , Ar o , Av o , Ay p ) (2.1.5) 
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The advantage of using equation (2.1.5) is that Ar p is almost independent of 
the semi-major axis a. In the range of a-values of interest, we may write 

Ar p = Ar p ( e > d > Ar o’ Av o> Ay o) (2.1.6) 

In this form Ar p depends only upon two parameters, the eccentricity e of 
the actual orbit and the true anomaly 6 at insertion. This greatly simplifies 
the statistical analysis of Ar p . 

2.2 Use of the Cumulative Distribution Function F 

An example of a cumulative distribution function 

F(Ar p ) = P(AR p <Ar p ) 

where AR = the random variable and A r = a particular value of AR , is shown 

P P P 

in Figure 2a. (For derivation of the distribution function see Chapter 2.3.) The 
use of this function can be demonstrated best by an example. 

Assume that the calculated perigee height (r Pcal - r e ) is 79 n.mi. (146 km). 
The allowed minimum perigee height (r - r e ) is 75 n.mi. (139 km), and it 
appears that a Go decision could be made. 1 ” However, the condition for making a 
correct Go decision is 


r ~ r >75 n.mi. 

p e — 


( 2 . 2 . 1 ) 


Using (2.1.4) and substituting AR p for Ar p (since we wish to make a state- 
ment of probability we will use the random variable AR p ) we can write 


Therefore, 


or 


(r 


Peal 


- r ) - AR > (r - r ) 

e ' n — v o : _ & ' 


79 n.mi. - AR > 75 n.mi. 

p — 

AR^ < 4 n.mi. (7.4 km) . 


( 2 . 2 . 2 ) 


From the cumulative distribution function in Figure 2b, we see that the 
probability for satisfying (2.2.2), i.e., making a correct Go decision, is: 

P(AR p 1 4 n.mi.) = F = 93%. (2.2.3) 


4 


F (^ r p) F ( dr ) F (dr 



0 5 n . mi . 


Figure 2a. Cumulative Distribution 


Function F of the Perigee Error 

AR F = F (Ar ) = P (AR < Ar ). 

K p P ~ p 


Figure 2b. With a Padding of 
4 Nautical Miles, the Probability 
for a Correct GO Decision i s 
93%. 


Figure 2c. A Padding of 7.4 
Nautical Miles Is Required for 
a 99.5% Probability of a Correct 
GO Decision. 


(km) 
d r p 
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Using the same function (see Figure 2c), for a 99.5% probability, the calcu- 
lated perigee height ( r - r ) has to be 

Peal e 

(r - r ) - 75 n.mi. > 7.4 n.mi. 

Peal e 


or 


(r 


Peal 


r e ) — 82.4 n.mi. 


(2.2.4) 


In other words, a padding of 7.4 n.mi. (13.7 km) is required to insure with a 
99.5% probability that the true or actual perigee height exceeds 75 n.mi. (139 km). 

The cumulative distribution function F used in the above example is valid 
f or e = 0.001, 6 = 180°, and no correlation between the random variables AR q , 
AV o , and Ar < 1 ) . The required padding for other insertion parameters and 
different correlations between AR q , AV q , and AT^ is analyzed in Chapter 3. 

2.3 Computation of the Cumulative Distribution Function F. 

Normally, the first step in an error analysis is to develop a variational 
equation expressing Ar p as a function of Ar o , Av o , and A 7q . This is done in 
Appendix B (Equation B.79), and it is shown that the variational equation is a 
second order equation in A r p of the form: 

a. Ar + a, Ar Ar + a Av Ar + a (Ar ) 2 

lp 2 o p 3op 4 v p 7 


^l Ar o + ^2 Av o + ' A 3 A ^o + A» Ar o Av o + / 3 S Av o A ^o 

(2.3.1) 


+ /3 6 Ar o A 7o + /3 ? (Ar o ) 2 + /3 g (Av q ) 2 + (AyJ 2 


+ terms of higher order 


where the coefficients a. (i =l,2,3,4)and /3. (i = 1,2, ..., 9) are functions of 
the insertion parameters r o , v o , and y o . 


^^Here AR q , AV q and AT o are random variables representing the errors in insertion radius, 

velocity magnitude, and flight path angle, to be distinguished from Ar q , Av q , and A which 

are values which AR , AV . and AT can assume. 

o’ o 
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Equation (2.3.1) has simple solutions for two special cases. For a cir- 
cular orbit, the coefficients of the linear terms vanish and the result is Equation 
(B.82) of Appendix B: 


A r 

p 




2 



( r o A T Q ) 2 • 


(B.82) 


For large e -values the second order terms may be neglected, and the re- 
sult is a simple linear expression in Ar o , Av o , and Ay . In particular, if 
0.005 £ e £ 0.05, equation (9) from reference 5 is valid 


Ar p % 


(2 - cos (9)Ar o + 2(1 - cos 9) 



(r o sin£)Ay o . (2.3.2) 


The range of e-values of interest in this report is 0 £ e £ 0.01. In the 
regionO £ e £ 0.005, the variational equation (2.3.1) confuses rather than aids 
the analysis. Therefore, the basic definition, equation (2.1.4), has been used. 


A, 


r ( r 

p V O 


+ Ar 


+ Ai 


y* + - 


r p Oo’ 


v o ,y o ). (2.1.4) 


It should be noted that this latter expression is exact. The value of equations 
(B.82) and (2.3.2) lies in their simplicity. Consequently, in their region of valid- 
ity, they have been used as a check against the numerical computations performed 
using (2.1.4). 

If AR q , AV q , and AT o are assumed to be normally distributed random vari- 
ables, then the cumulative distribution function of perigee error 

AR = r (r + AR , V + AV , y + AF ) - r ( r , v , y ) (2.3.3) 

p p v o o o o o o 7 p v o o / o / ' ' 


can be written in integral form. However, this integral cannot be expressed in 
terms of known or tabulated functions, but must be evaluated by numercial tech- 
niques - either a Monte Carlo approach or numerical integration. The latter 
approach has been used in this analysis, and is described in Appendix B. 

If the error in perigee is expressed in terms of the orbital elements 
a , e, d 


AR 


= AR p (a, 


e, 9, AR , AV , Ar ) 

o’ o’ n / 


(2.3.4) 


( 2 ) 


Again, since we are dealing with probability distributions we will use random variable 


notations. 
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-20 -15 -10 -5 0 5 10 15 20 25 (km) 

I i l I I l I l i I i I I 1 I I 1 I I 1 I I 1 1 I I Arp 

- 10 -5 0 5 ,0 (n. mi .) 

Figure 3. The Cumulative Distribution Function F (Ar p ) Varies Insignificantly if the Semi-Major 
Axis a Is Varied ±54 Nautical Miles (+100 km) Which Is Approximately Six Times the Variation 
in a for this GO, NO-GO Analysis. Insertion Parameters: Eccentricity e = 0.01, True Anomaly 

at Insertion 6 = 225°. 


we see in Figure 3 that AR p is almost independent of the semi -major axis a for 
a-values of interest in the Go, No-Go analysis. Figure 3 shows the cumulative 
distribution functions for AR p for the following values of a: 

3490 n.mi., 3544 n.mi., and 3598 n.mi. 

(6463 km, 6563 km, 6663 km) 

Although a varies ± 54 n.mi. (±100 km), or approximately six times the 
expected variation in a for this analysis, the cumulative distribution functions 
vary insignificantly. We may therefore neglect the variation in a , and use an 
average value of 3544 n.mi. (6563 km). This will greatly simplify the analysis. 
Hence, 


AR = AR (e, 0, AR , AV , Ar ). (2.3.5) 

p p v ’’ooo / 

A computed probability density function of perigee error, f (Ar ) for an 
actual circular parking orbit (e = 0) for both correlated and uncorrelated inser- 
tion errors AR q , AV q , and AP o is shown in Figure 4a. The computed cumulative 
distribution function 
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F (Ar 



Figure 4a. Probability Density Function f of Perigee Error AR for a Circular Orbit. Insertion 

P 

Parameters: Semi-Major Axis a = 3544 n.mi., Eccentricity e - 0. 



Figure 4b. Cumulative Distribution Function F for the Above Probability Density Function f. 








(2.3.6) 


F 


F( A r p ) 


• A 'P 

f (x) dx 


is shown for the same conditions in Figure 4b. It is of interest to note that the 
perigee error does not have a mean of zero and is not normally distributed. 

Figures 5 through 8 show cumulative distribution functions for eccentrici- 
ties of e = 0.0005, e = 0.001, e = 0.005, and e = 0.01. In each instance there are 
curves for both correlated and uncorrelated insertion errors. For uncorrelated 
errors, curves are plotted for true anomalies at insertion of 9 = 0°, 90°, and 
180°. The F -curves are even functions of 6 for uncorrelated errors and hence 
are also valid for 9 = -90°. For correlated errors, the curves are no longer 
even functions of 9, and have therefore been plotted for 9 = 0 ° and 9 = 180° ±45° 
In all of these figures, the "3 sigma" limits for the insertion errors have been 
taken as 


3o^ r - 2.4 n.mi. (4.44 km) 

o 

3ct^ v = 16 ft/sec (4. 87 tn/sec) 

3c7 Ay 0 = 0.16° (2.79 m rad) . 


2.4 Verification of the Numerical Integration 

It has already been mentioned that in its region of validity equation (2.3.2) 
is a very good approximation for the perigee error. This can be seen in Figures 
9a and 9b where the cumulative normal distribution function using (2.3.2) is 
compared with the cumulative distribution function of perigee error obtained 
by numerical integration. In both figures, the semi-major axis was taken to be 
3544 n.mi. (6563 km) and the eccentricity of the parking orbit e = 0.01. In Fig- 
ure 9a, 6 = 225° and the coefficients of correlation of the errors at insertion 
AR q , AV o , and AF q are all equal to 0.75. It can be seen that the agreement be- 
tween the two methods is very good. In Figure 9b, 9 = 180°, and the errors 
are uncorrelated. Again the agreement between the two methods is very good. 

In Appendix B, Figures 19 and 20, the required padding for 99.5% probabil 
ity is compared using both numerical integration and equation (2.3.2). Again, 
there is good agreement. 
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Figure 5a. Cumulative Distribution Function F of Perigee Error AR . Insertion Parameters 
a - 3544 n.mi. (6563 km), e - 0.0005, True Anomaly at Insertion 0 - 0°, 135°, 225°. Inser- 
tion Errors Are Positively Correlated, Coefficients of Correlation ~ +0.75. 



Figure 5b. Cumulative Distribution Function F of Perigee Error AR . Insertion Parameters 
□ - 3544 n.mi. (6563 km), e = 0.0005, 0 - 0°, 90°, 180°. Insertion Errors Are Uncorrelated 




F (Ar p ) n ° T1 F ( Ar p ) 



igure 6a. Cumulative Distribution Function F of Perigee Error AR^. Insertion Parameters: 

135°, 225°. Insertion Errors Are Positively 


- 3544 n.mi. (6563 km), e - 0.001, 8 - 0 C 
orrelated, Coefficients of Correlation ~ +0.75. 



Figure 6b. Cumulative Distribution Function F of Perigee Error AR^. Insertion Parameters: 
a = 3544 n.mi. (6563 km), e = 0.001, 8 - 0°, 90°, 180°. Insertion Errors Are Uncorrelated. 
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uo 


25 (km) 


Figure 7a. Cumulative Distribution Function F of Perigee Error AR . Insertion Parameters: 
a - 3544 n.mi. (6563 km), e 3 0.005, 6 - 0°, 135°, 225°. Insertion Errors Are Positively 
Correlated, Coefficients of Correlation - +0.75. 


Figure 7b. Cumulative Distribution Function F of Perigee Error AR . Insertion Parameters: 


a H 3544 n.mi. (6563 km), e ~ 0.005, 6 ~ 0°, 90°, 180°. Insertion Errors Are Uncorrelated 












Figure 8a. Cumulative Distribution Function F of Perigee Error AR^. Insertion Parameters: 
a = 3544 n.mi. (6563 km), e = 0.01, 6 = 0°, 135°, 225°. Insertion Errors Are Positively 
Correlated, Coefficients of Correlation “ +0.75. 


LEGEND 


Figure 8b. Cumulative Distribution Function F of Perigee Error AR p . Insertion Parameters 
a = 3544 n.mi. (6563 km), e “ 0.01, 9 - 0°, 90°, 180°. Insertion Errors Are Uncorrelated 








Figure 9a. The Computed Cumulative Distribution Function F of Perigee Error AR^ Agrees with 
the Cumulative Normal Distribution Approximation for Large Eccentricities. Insertion Parame- 
ters: a - 3544 n.mi. (6563 km), e = 0.01, 9 - 225°. Insertion Errors Are Positively Corre- 
lated, Coefficients of Correlation = +0.75. 



Figure 9b. The Computed Cumulative Distribution Function F of Perigee Error AR p Agrees With 
the Cumulative Normal Distribution Approximation for Large Eccentricities. Insertion Parame- 
ters: a - 3544 n.mi. (6563 km), e - 0.01, 6 - 180°. Insertion Errors Are Uncorrelated. 


15 




3. THE GO, NO-GO DECISION 


An incorrect Go decision (one where the true or actual perigee height is less 
than 75 n.mi.) could endanger the lives of the astronauts. Therefore, a 99.5% 
probability is required for a correct Go decision (reference 1). Since there will 
always be an error Ar p in the calculation of the perigee height, a padding must 
be added to the minimum required perigee height of 75 n.mi. (139 km) in order 
to have a 99.5% probability that the actual perigee height exceeds 75 n.mi. The 
required padding must be equal to or greater than that value of A r p for which 
F(Ar p ) = 0.995. Hence, it can be read off the cumulative distribution curves in 
Chapter 2 at the point where F = 0.995. 

In this chapter, the required padding is analyzed and is seen to depend sig- 
nificantly on the eccentricity e, the true anomaly at insertion 9 , and the corre- 
lation coefficients between the insertion errors. As discussed in paragraph 2.1, 
the required padding depends insignificantly on the semi-major axis a . Since it 
is practically impossible to determine the actual coefficients of correlation 
between the insertion errors, the calculations have been carried out for coeffi- 
cients of correlation of 0, ±0.75, and ±0.9. That is, all possible combinations of 
signs have been considered with the exception of those for which the covariance 
matrix becomes singular. Throughout the analysis, it has been assumed that the 
insertion errors AR o , AV o , and AF q are normally distributed random variables 
having means of zero and "3 sigma" values of: 

3o-^ r = 2.4 n.mi. (4.44 km) 

3<t^ v - 16 ft/sec (4.87 m/sec) 

3cr Aro = 0. 16° (2.79 m rad) 

Table 1 shows the maximum required padding for eccentricities from 0 
through 0.01. Also shown are the associated true anomalies at insertion and the 
coefficients of correlation for the insertion errors. 

Figures 10 through 14 show the required padding for a 99.5% probability of 
a correct Go decision. The figures are for eccentricities of 0, 0.0005, 0.001, 
0.005, and 0.01, and the curves are plotted versus true anomaly at insertion 6 . 

In each figure there are two sets of curves — those for correlation coefficients 
of 0.75 and those for correlation coefficients of 0.9. 

It is of interest to note that for eccentricities in the range e < 0.001, the 
maximum required padding occurs when the insertion errors are uncorrelated. 
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Table 1 


Maximum Required Padding for a 99.5% 
Probability for a Correct Go-Decision 


Eccentricity 

Max. Required 
Padding 

True anomaly 
at insertion 
(degrees) 

Coefficients of Correlation 

^Ar Av 

o o 

^Av Ay 

o ' o 

0 

0 

<5 

Q. 

(n.mi.) 

(km) 

0 

n 


All values 

0 

0 

0 


4.8 

8.9 

180 


0 

0 


7.5 

13.8 

180 

0 

0 

0 

0.005 

17.1 

31.7 

135 

+0.9 

-0.9 

-0.9 

225 

+0.9 

+0.9 

+0.9 

0.01 

17.2 

31.8 

135 



-0.9 

225 

+0.9 

+0.9 

+0.9 


For eccentricities in the range e > 0.005, high correlation gives the worst case. 
The maximum required padding does not increase substantially if e is increased 
above 0.005. This result obtained from numerical integration also agrees with 
equation (2.3.2) which is valid in the region 0.005 < e < 0.05, and is independent 
of e. 

4. CONCLUSIONS 

4.1 The Optimum Orbit 


From the standpoint of the Go, No-Go decision, a circular parking orbit is 
most desirable for two reasons: 

1. For a given insertion energy, the perigee height is the greatest, giving 
the largest margin between the actual perigee height and the required 
minimum. 
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Figure 10a. Required Padding as a Function of True Anomaly at Insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r pmjn - r e ) of 
75 n.mi. (139 km). Insertion Parameters: a = 3544 n. mi . (6563 km), e = 0 (Circular Orbit). 

Correlated (All Possible Combinations of Signed Values of 0.75) and Uncorrelated Insertion 
Errors. 



Figure 10b. Required Padding as a Function of True Anomaly at Insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum ( r Pmir| ” r e ^ 

75 n.mi. (139 km). Insertion Parameters: a = 3544 n.mi. (6563 km), e - 0 (Circular Orbit). 

Correlated (All Possible Combinations of Signed Values of 0.90) and Uncorrelated Insertion 
Errors. 






J TRUE ANOMALY 
AT INSERTION 0 
360 (degrees) 


Figure 11a. Required Padding as a Function of True Anomaly at Insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r ~ r e ) of 
75 n.mi. (139 km). Insertion Parameters: a “ 3544 n.mi. (6563 km), e = 0.0005. Correlated 
(All Possible Combinations of Signed Values of 0.75) and Uncorrelated Insertion Errors. 



— TRUE ANOMALY 
— 1 AT INSERTION 0 
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Figure lib. Required Padding as a Function of T rue Anomaly at I nsertion in Order to I nsure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r “ r ) of 

Pmin e 

75 n.mi. (139 km). Insertion Parameters: a - 3544 n.mi. (6563 km), e = 0.0005. Correlated 
(All Possible Combinations of Signed Values of 0.90) and Uncorrelated Insertion Errors. 





Figure 12a. Required Padding as a Function of True Anomaly at Insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r “ r ) of 

P m i n ” 


75 n.mi. (139 km). Insertion Parameters: a - 3544 n.mi. (6563 km), e - 0.001. Correlated (All 


Possible Combinations of Signed Values of 0.75) and Uncorrelated Insertion Errors. 



TRUE ANOMALY 
AT INSERTION 9 
(degrees) 


Figure 12b. Required Padding as a Function of True Anomaly at Insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r p “ r e ) of 
75 n.mi. (139 km). Insertion Parameters: a - 3544 n.mi. (6563 km), e ” 0.001. Correlated (All 
Possible Combinations of Signed Values of 0.90) and Uncorrelated Insertion Errors. 
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Figure 13a. Required Padding as a Function of True Anomaly at insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r “ r^) of 
75 n.mi . (139 km). Insertion Parameters: a = 3544 n.mi. (6563 km), e = 0.005. Correlated (All 
Possible Combinations of Signed Values of 0.75) and Uncorrelated Insertion Errors. 



Figure 13b. Required Padding as a Function of True Anomaly at Insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r “ r # ) of 
75 n.mi. (139 km). Insertion Parameters: a - 3544 n.mi. (6563 km), e — 0.005. Correlated (All 
Possible Combinations of Signed Values of 0.90) and Uncorrelated Insertion Errors. 




Figure 14a. Required Padding as a Function of True Anomaly at Insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r ” r ) of 
75 n.mi. (139 km). Insertion Parameters: a - 3544 n.mi. (6563 km), e - 0.01. Correlated (All 
Possible Combinations of Signed Values of 0.75) and Uncorrelated Insertion Errors. 



Figure 14b. Required Padding as a Function of True Anomaly at Insertion in Order to Insure with 
a 99.5% Probability that the Actual Perigee Height Exceeds the Allowed Minimum (r “ r e ) of 
75 n.mi. (139 km). Insertion Parameters: a - 3544 n.mi. (6563 km), e - 0.01. Correlated (All 
Possible Combinations of Signed Values of 0.90) and Uncorrelated Insertion Errors. 
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2. As shown in Chapter 3, the perigee error is smallest. 

Since the circular orbit cannot be achieved in actual flight, the optimum or- 
bit is therefore one which has the smallest eccentricity. In addition, as shown in 
Chapters 2 and 3, insertion at a true anomaly of 9 - 0° is desirable in order to 
minimize the error in perigee Ar p . 

4.2 The Required Padding 


In Chapter 3 it was shown that the required padding depends significantly on 
the eccentricity e of the parking orbit, increasing with e in the region 0 < e < 
0.005, and remaining nearly constant in the region 0.005 < e < 0.01 (Figures 10 
through 14). In order to use the smaller required padding corresponding to 
smaller e-values, the insertion ship must have the capability of determining the 
eccentricity of the orbit with sufficient accuracy. This can be clarified with an 
example. 

The maximum required padding for an eccentricity of e = 0.001 is 7.5 n.mi. 
(14 km) according to Figures 12a and 12b in Chapter 3. In order to use this 
value, however, the ship must be able to verify that e < 0.001 for the actual or- 
bit. The calculated eccentricity e cal is a function of the measurement errors. 
Figures 15a and 15b show the cumulative distribution function H (Ae) of the error 
in eccentricity AE for orbits with 


e 

0.001 



3cr = 

r o 

2. 4 n.mi. 

(4.44 

km) 

3 a = 

v o 

16 ft/sec 

(4.87 

m/sec) 

3c Vo = 

0. 16° (2. 

79 mrad) 


For the derivation and calculation of H(Ae) see Appendix C. 

It is seen from these curves that the insertion ship cannot verify that 
e < 0.001. Therefore, the small padding of 7.5 n.mi. (14 km) cannot be used. 

The cumulative distribution function H(Ae) for an actual orbit with e = 0.005 
is shown in Figures 16a and 16b. It is evident from these figures that even an 
eccentricity as large as e = 0.005 cannot be determined accurately enough by 
the tracking ship. Therefore, the smaller required padding corresponding to 
the lower e -values cannot be used. 
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Figure 15a. Cumulative Distribution Functions of the Error in Eccentricity for Various Values of 
True Anomaly at Insertion. Insertion Parameters: a - 3544 n.mi. (6563 km), e - 0.001. Un- 
correlated Insertion Errors. 



Ae 


Figure 15b. Cumulative Distribution Functions of the Error in Eccentricity for Various Values 
of True Anomaly at Insertion. Insertion Parameters: a ~ 3544 n.mi. (6563 km), e ” 0.001. 

Correlated Insertion Errors. (p Ar . - 0.9, P Av Ay ~ ^ Ar Ay “ ±0-9) 





H (Ae) 


1 .0 


0 . 01 = 

- 0.002 



Figure 16a. Cumulative Distribution Functions of the Error in Eccentricity for Various Values 
of True Anomaly at Insertion. Insertion Parameters: a - 3544 n.mi. (6563 km), e “ 0.005. 

Uncorrelated Insertion Errors. 


LEGEND 


Figure 16b. Cumulative Distribution Functions of the Error in Eccentricity for Various Values 
of True Anomaly at Insertion. Insertion Parameters: a - 3544 n.mi. (6563 km), e - 0.005. 

Correlated Insertion Errors (p/± f = +0.9, ^ y ~ p^ r {± y = ±0.9) 






As previously mentioned, the required padding does not increase significantly 
for e > 0.005. It is therefore recommended that the maximum padding of 17 n.mi. 
(32 km) for e > 0.005 be used as a required padding. This value insures with a 
99.5% probability that the actual perigee height exceeds the required minimum 
perigee height for the following cases: 

1. All eccentricities in the region 0 < e < 0.01. 

2. All values of true anomaly at insertion. 

3. All values of correlation coefficients p between the insertion errors in 
the region 0 < |p| < 0.9. 

4.3 Recommendations 

The results of this analysis may be briefly stated in the form of the follow- 
ing recommendations: 

1. A padding (calculated perigee height > required minimum perigee 
height + padding) of 17 n.mi. (32 km) should be used in order to insure 
with a 99.5% probability that a positive Go decision is correct. 

2. The near-earth Apollo parking orbit should have the smallest possible 
eccentricity in order to achieve the largest possible margin between 
the actual perigee height and the required minimum perigee height. 

3. Insertion of the Apollo spacecraft should occur at a true anomaly as 
close as possible to 6 = 0° in order to minimize the perigee height 
error. 
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GLOSSARY AND DEFINITION OF SYMBOLS 


Random vector . Let Vj, V 2 V be p random variables. Then the (pxl) 

vector V, 



is a random vector. 

Multivariate Normal Distribution. Let the p -dimensional random vector V have 
the probability density function 


|nr* ( 2 tt )~'2 e -* (V-^) T n- 1 (V-n) 

where 7? is the vector of constants 



and n is a (p x p) positive definite matrix. Then V has a nonsingular multi- 
variate normal distribution with mean vector £(V) = 7 j and covariance matrix 

s[(v - ij)(v- 7}) t ] =n. 

Uncorrelated The random variables V x and V 2 are said to be uncorre- 
lated if the coefficient of correlation 0 „ between them is 

v 1 V 2 


zero. 


Error 


The measured or calculated value of a quantity minus its 
actual or true value 


Semi -major axis of the parking orbit 
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Element in the i th row and j th column of the orthogonal 
matrix B relating the normal vector Y (components are 
uncorrelated) to the normal vector X (components are 
correlated) such that Y = BX 

The (3 x 3) diagonal matrix relating the normal vector W to 
the normal vector X such that W = AX. The matrix A has 
the diagonal elements <r A , a A , and <x A . 

The (3 x 3) orthogonal matrix relating the normal vector 

Y (components are uncorrelated) to the normal vector X 
(components are correlated) such that Y = BX. The matrix 
B is also the transpose of the matrix C . 

Element in the i th row and j th column of the (3 x 3) ortho- 
gonal matrix C relating the normal vector Y (components 
are uncorrelated) to the normal vector X (components are 
correlated) such that X = CY. The matrix C is also the 
transpose of the matrix B . 

The (3 x 3) orthogonal matrix relating the normal vector 

Y (components are uncorrelated) to the normal vector X 
(components are correlated) such that X = CY. The matrix 
C is also the transpose of the matrix B. 

An orthogonal matrix (3 x 3) which diagonalizes the corre- 
lation matrix P lf i.e., C i T P jCj is a diagonal matrix. 

An orthogonal (3 x 3) matrix which diagonalizes the corre- 
lation matrix P 2 , i.e., C 2 T P 2 C 2 is a diagonal matrix. 

An orthogonal (3 x 3) matrix which diagonalizes the corre- 
lation matrix P 3 , i.e., C 3 T P 3 C 3 is a diagonal matrix. 

A (3 x 3) diagonal matrix having \ lf \ 2 , and X. 3 as diagonal 
elements. D = E T PE where E is an orthogonal matrix and 
P is the covariance matrix of X (correlated components). 

A calculated value of eccentricity. 

e . = e( r + Ar , v + Av , y + Ay ) 

The actual or true value of eccentricity. 


e(Xj) 

E 

f 2 (y 2 ) 

f 3 (y 3 ) 

f (yi* y 2 * y 3 

f 

Fi(Ar o ) 

F 

g( x ) 

g(x,, x 2 , x 3 

G(x) 

H (A e) 

I 

J (Xj, X j * X 3 


The expected value of the random variable X x . Sis the 
expectation operator. 

A (3 x 3) orthogonal matrix which diagonalizes the corre- 
lation matrix P, i.e., D = E T PE is a diagonal matrix having 
k v \ 2 , and k 3 as diagonal elements. 

Probability density function of the normally distributed 
random variable Y . 

Probability density function of the normally distributed 
random variable Y 2 . 

Probability density function of the normally distributed ran- 
dom variable Y 3 . 

Joint probability density function of Yj, Y 2 , and Y 3 where 
Yj, Y 2 , and Y 3 are independent normal random variables. 

f (yi»y 2 .y 3 ) = f i(yi) f 2 (y 2 ) f 3 (y 3 )- 

The probability density function of the error in perigee AR . 
We have f = F'(Ar p ) = f(Ar ) ' 

Cumulative distribution function of AR 

O 

The cumulative distribution function of the error in perigee 
A R p where F = P(AR p < Ar p ) = F(Ar p ) 

The probability density function of the standardized normal 
random variable X . 

The joint probability density function of X 3 , X 2 , and X 3 
where Xj, X 2 , and X 3 are components of the normal vector 

X. 

The cumulative distribution function for the standardized 
normal random variable X . 

The cumulative distribution function for the error in ec- 
centricity AE. 

The (3 x 3) identity matrix. 

The Jacobian of the transformation between Xj , X 2 , X 3 and 
Yi, Y 2 ,Y 3 . 
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K(A 6) 

1 


L 


m 


mod 


( 2N + 1) 


Pi, 


J (Xj, X 2 , X 3 ) - mod 


3 v« 

- i 

3v. 

^y 

B x j 

3 x 2 

B x 

3y 2 

Sy 2 


d Xj 

3 x 2 

B x 

^y 3 

3 y 3 

By 

B x x 

B x 2 

B x 

absolute value of the 

expression 


following it. 

The cumulative distribution function for the error in true 
anomaly A0. 


A constant given by (see m also). 
/2 r \ 


1 = 


°7W„ + 2cr Ar f 


(m > 1) 


Upper and lower limit (in number of standard deviations) 
for the approximation to the range of each independent 
normal random variable Yj (i = 1,2,3). 

A constant given by (see 1 also) 



where (m > 1) 

An operator which means the absolute value of the expression 
following it. 

Total number of subdivisions (extending from -Lo- to +L cr ) 
for the approximation to the range of each independent 
normally distributed random variable Y A (i = 1, 2, 3). 

The probability that the discrete approximation Yd to the 
normally distributed random variable Y £ assumes the value 
y! where y ' is the midpoint of the interval 
(y i k - Ay j , yj k + Ay 4 ) . We have 
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P(y'„ -ay, < V, < y\ u *iy,) a p(V, = y :j = Pi> 


where 


(i = 1, 2, 3) and (k = 1, 2, 3, .... (2N + 1)) 


Also 


q i, 

P ik “ ~ ( see q ik s i> 


P 


jki 


The joint probability that the discrete random variable ap- 
proximations Yj, Y 2 , Yg to the independent normal random 
variables Y lt Y 2 , Y 3 assume the values , y' 2 , y 3 . We 
have J 


JH 


p ( Y i = 


Y' = 
1 2 


y 2l 


Y' = 

1 3 


y 3.) 




where (j, k, 1 = 1, 2, 3, . . . , 2N + 1) 

p(Yj “ y- ) The probability that Y! (i = 1, 2, 3) takes the value y! 

k/ (k = 1, 2, 3, . . 2N + 1). P denotes "probability" here. 

P Correlation matrix or covariance matrix of the normal vec- 

tor X when p = p =p = p, (0 < p < 1). 

’ X j X 2 ^ X 2 X 3 ^ x l x 3 r — r 

A real symmetric positive definite square matrix of order 3. 

Pj Correlation matrix or covariance matrix of the normal vec- 

tor X when p xjx2 = -p, p xjx? = -p, p X2X3 = p (0 < p < 1). 
A real symmetric positive definite square matrix of order 3. 

P 2 Correlation matrix or covariance matrix of the normal vec- 

tor X when p xix2 = -p, p XlX 3 = p, p X2X3 = -p, (0 < p < 1). 
A real symmetric positive definite square matrix of order 3. 

P 3 Correlation matrix or covariance matrix of the normal vec- 
tor X when = p, p Xl x 3 = ~P> P* 2 * 3 = "P. (0 < p < 1). 

A real symmetric positive definite square matrix of order 3. 

q A Value of the normal probability density function of Y t (i =1, 

2, 3) at the point y' (k = 1, 2, 3, . . 2N + 1), i.e., 

1 k 
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0 


r 


e 


r 

p 


r 

Pmin 


r 

o 


s. 


v 

o 


w 


x- space 


x 


l 


X 


2 


(- 1 


yi 






2 /2a y \ 


(see also p. ). 

1 k 

The covariance matrix of the normal vector Y . We have 
Q = S(YY t ) = BPB T = C T PC. 

Radius of the earth. 

Perigee radius. 

Calculated perigee radius (from the ship's tracking data). 
Required minimum perigee radius. 

Magnitude of the radius vector at insertion. 

A normalizing sum, i.e., 

2N+ 1 

s > = 23“ « (i = *' 2 ' 3> 

k= 1 

(see also p t ). 

The magnitude of the velocity vector of insertion. 

Normal random vector having the components AR o , AV o , and 

Ar . 

o 

Refers to the correlated "space" of the normal random vari- 
ables Xj, X 2 , and X 3 . 

Normalized eigenvector corresponding to the eigenvalue 

\ 1 = 1 + 2p (0 < p < 1). 

Normalized eigenvector corresponding to the eigenvalue 
X 2 = 1 ~ P (0 < p < 1). 

The dot product of the normalized eigenvector with the 
normalized eigenvector x 2 . 
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The normalized eigenvector corresponding to the eigenvalue 
^3 = 1 " P (0 < p < 1). 

Normal random vector having correlated standardized normal 
random variables X x , X 2 , and X 3 as components. 

Standardized normal random variable. A component of X . 

Standardized normal random variable. A component of X. 

Standardized normal random variable. A component of X. 

Discrete random variable approximations to the standardized 
normal random variables X x , X 2 , and X 3 . 

A standardized normal random variable. 

Refers to the uncorrelated "space" of the uncorrelated nor- 
mal random variables Y v Y 2 , and Y 3 . 

Normal random vector having the uncorrelated normal random 
variables Y x , Y 2 , and Y 3 as components. 

Values that the discrete random variable approximations 
Y x , Y 2 , Y 3 to the uncorrelated normal random variables 
Yj, Y 2 , Y 3 can take (j, k, 1 = 1, 2, 3, . . ., 2N + 1). 

Uncorrelated normal random variables. Components of the 
normal vector Y . 

Discrete random variable approximations to the uncorrelated 
normal random variables Y v Y 2 , and Y 3 . 

Normal random vector having the uncorrelated standardized 
normal random variables Z v Z 2 , and Z 3 as components. 

Uncorrelated standardized normal random variables. Com- 
ponents of the normal random vector Z. 

Coefficients in variational equation for A r p . Functions of 

r o» v o> 7 0 - (i = 2 » 3 > 4 ) 

Coefficients in variational equation for &r p . Functions of 
r o , v o , and y Q . (i = 1, 2, . . ., 9). 


A e A value for the error in eccentricity AE due to metric track- 

ing errors by the insertion ship. We have Ae = e cg . - e . 

AE A random variable representing the error in eccentricity. 

AE has the cumulative distribution function H(Ae). 

A y A value for the error in flight path angle at insertion AT o 

where AT q is normally distributed. 

AT o A normal random variable representing the error in flight 

path angle at insertion. 

AP A discrete random variable approximation to the error in 

flight path angle at insertion Ar o (Ar o is normally distributed). 

Ar p A value for the error in perigee radius AR p . 

AR p A random variable representing the error in perigee radius. 

AR p has the cumulative distribution function F(Ar p ), and 
probability density function f (Ar p ). 

AR p A discrete random variable approximation to the error in 

perigee radius AR p . 

Ar A value for the error in insertion radius AR where AR is 

o o o 

normally distributed. 

A<9 The error in true anomaly at insertion. A value A© can as- 

sume. 

A© A random variable representing the error in true anomaly of 

the spacecraft at insertion. 

AR o A normal random variable representing the error in insertion 

radius. 

AR° A discrete random variable approximation to the normal 

random variable AR o where AR o represents the error in 
radius at insertion. 

Av o A value for the error in the magnitude of the insertion velocity 

AV o which is a normal random variable. 
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AV 0 

av; 


Ayi 


y 

' o 


e 


k 


^l» "k 2 > 


A 


P 


^Ar 0 Av 0 


^Av 0 A 7o 


^Aro Ay c 


A normal random variable representing the error in the 
magnitude of the velocity at insertion. 

A discrete random variable approximation to the normal 
random variable AV o which represents the error in inser- 
tion velocity magnitude. 

A subdivision of the range of the normal random variable 
approximation for Y . (i = 1, 2, 3). We have 

Ay. = (— 7 ) 0 - where a is the standard deviation of Y. . 

Flight path angle at insertion measured from the normal to 
the radius vector at insertion, pointing in flight direction, 
against the direction of motion (i.e., clockwise if motion is 
counterclockwise), such that the direction at y Q = 90° is 
radially away from the center of attraction. 

True anomaly of the spacecraft at insertion into parking orbit 
the angle measured from the perigee radius to the radius vec- 
tor at insertion in the direction of motion. 

A variable in the secular (characteristic) equation of P . 

Eigenvalues of the covariance matrix P or elements of the 
diagonal matrix C T PC, the covariance matrix of Y . We have 
= * y \ * ^2 = CT y 2 2 . ^ A 3 = a y 2 . 

Covariance matrix of the normal random vector W having 
correlated normal random variables AR o , AV 0 , and AT o as 
components. 

Gravitational constant of the earth (398603.2 km 3 /sec 2 ). 
Coefficient of correlation ( \p\ ± 1). 

Coefficient of correlation between AR o and AV o . 

Coefficient of correlation between AV o and Ar o . 

Coefficient of correlation between AR 0 and AP 0 . 
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p 


x 1 x 2 


Coefficient of correlation between and X 2 
are standardized normal random variables) . 


(X x and X 2 


P X 2 x 3 Coefficient of correlation between X 2 and X 3 . (X 2 and X 3 

are standardized normal random variables). 

P Xi x 3 Coefficient of correlation between X x and X 3 .(X x and X 3 

are standardized normal random variables) . 


cr , o- , cr Standard deviations of the uncorrelated normal random vari- 

y i y 2 y 3 

ables Y x , Y 2 , and Y 3 . 

crj 2 , cr y 2 , p - y 2 Variances of the uncorrelated normal random variables Y x , 
Y a , Y 3 . Diagonal elements of C T PC. 


a 


Standard deviations of the correlated normal random variables 

( r o + AR o)» < V o + AV o)» 311(1 (y o + Ar o>- 


°"Ar > °Av „ > a S.y 


Standard deviations of the correlated normal random variables 
AR o , AV o , and Ar o . 


cr ^ 2 , cr A 2 , cr£ y Variances of the correlated normal random variables AR o , 

AV o , and Af o . Diagonal elements of A, the covariance matrix 
of W . 


^Ar- Av, 


^Avo A y 0 


Cr \r 0 bry 0 


<7 


X 1 x 2 


Covariance between the normal random variables AR„ and AV . 
Off diagonal element of A , the covariance matrix of W . We 

haVe CT Ar 0 Av 0 = ^AroAv^Ar^Av, • 

Covariance between tlje normal random variables AV 0 and 
Ar o . Off diagonal element of A , the covariance matrix of W . 
Wehave o- Avo a 7o = PAv 0 Ay 0 CT Av 0 ^A 7o • 

Covariance between the normal random variables AR and 

O 

Ar o . Off diagonal element of A, the covariance matrix of W. 
Wehave a AroAyo = P A r 0 A 7o CT Ar 0 CT Ay 0 • 

Covariance between the standardized normal random vari- 
ables X l and X 2 . Off diagonal element of P , the covariance 
matrix of X (also called the correlation matrix). We have 
ct „ = p v cr cr = p v „ since X. and X, are 

standardized variables. 
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Covariance between the standardized normal random variables 
X 2 and X 3 . Off diagonal element of P, the covariance matrix 
of X . We have o- = p a a = p since X, and 

x 2 x 3 r *2 x 3 x 2 x 3 r *2*3 2 

X 3 are standardized variables. 

Covariance between the standardized normal random variables 
Xj and X 3 . Off diagonal element of P, the covariance matrix 
of X . We have a = p a cr = p since both X. 

and X 3 are standardized variables. 



APPENDIX A 


EXPRESSION FOR THE PERIGEE RADIUS IN TERMS OF THE INSERTION 
PARAMETERS 

The basic Keplerian two-body equations of motion which relate the insertion 
parameters r o , v q , and y o to the orbital elements a, e, and 9 (see Figure 1) 
are: 


(1 - e 2 ) 


■'G-t) 

(A.l) 

2 2 2 
v * r z cos y 

o o 7 o 

fJL a 

(A.2) 

a ( 1 - e 2 ) 
1 + e cos 6 

(A.3) 


where r , v o , y o are respectively the magnitude of the radius vector at insertion, 
the magnitude of the velocity vector at insertion, and the flight path angle at in- 
sertion. The semi-major axis is a, the eccentricity of the parking orbit is e, 
and the true anomaly of the spacecraft at insertion is 6. The gravitational con- 
stant iS (A. 

The perigee radius r p is the value of r o when 6 = 0°. From (A.3) we have: 

r p = a (1 - e) . (A.4) 


Rewriting (A.l) 


a 



(A.5) 


Substituting (A.5) into (A.2) and solving for e we have 


e 



y* 


+ 




cos 2 y 

' o 


(A.6) 
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where the positive sign is taken by definition. Using (A.5) and (A.6) we can write 
(A.4) as: 



APPENDIX B 


THE CUMULATIVE DISTRIBUTION FUNCTION OF PERIGEE ERROR 
BY NUMERICAL INTEGRATION 

The perigee radius r p can be expressed as a function of the insertion par- 
ameters r o , v o , and y o (Appendix A) 



(B.l) 


In calculating these insertion parameters from measurements, the errors 
Ar o , Av o , and Ay o are introduced. Therefore, the calculated perigee radius 
r will deviate from the actual perigee radius r by the amount A r . 

Peal fl £ P J p 


Peal 


r (r + Ar , v + Av , y + Ay ) - r (r , v , y ) 

p v o o’ o o’ / o / o / p v o’ o’ 'o' 


(B.2) 


Let us denote by AR q , AV q , and AF q random variables representing the er- 
rors in the calculation of the insertion parameters/ 1 } Further, let us assume 
that the random vector 


W = 



(B.3) 


has a trivariate normal distribution with zero mean vector and covariance 
matrix A 


^It should be noted that Ar is a value which the random variable AR can assume. The proba- 
bility that AR q is less than or equal to Ar Q is a function of Ar Q . Mathematically, we write 

F, (Ar ) = P (AR < Ar ) 

1 v o 7 v o — o 7 

where (Ar Q ) is the cumulative distribution function of AR q . 
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f < 

CT Ar 0 Av 0 

> q 

►n 

o 

* 

A = 

a Ar 0 Av 0 

1 

CT A 2 v q 

o 

> 

b 1 


£ 

n 

0 

$ 

o 

^ q 

< 

0 

% 

O 



J 


(B.4) 


where cr^ 2 ,a£ v , and cr£ y are the variances of AR q ,AV o , and Ar o respectively, 
and cr» * , cr\ * y , and cr. are the covariances between these random 

A* r o L ' v o o ' o r 0 l \ / 0 

variables . ( 2) 

Then W has the density function ( 3 } 


-A - — w T A' 1 w 

(277) 2 |A| 2 e 2 


(B.5) 


Since W has a trivariate normal distribution, AR , AV , and Ar are 

o O' o 

normally distributed random variables, and 


AR = r (r +AR,v +AV,y +Ar)-r (r , v , y ) (B.6) 

p p v o o’ o o’ / o o' p v o’ o’ 'o' ' ' 


is a random variable with a cumulative distribution function F (A r p ) and a prob- 
ability density function f (Ar p ). Although F(Ar p ) can be expressed in integral 
form, the integral cannot be expressed in terms of known or tabulated functions. 

It must be evaluated by numerical techniques — either a Monte Carlo approach 
or numerical integration. The latter approach will be described in the following 
paragraphs. In an effort to simplify notation we will use matrix algebra wherever 
possible. 

Since AR o , AV o , and AT o are jointly distributed with a trivariate normal dis- 
tribution, they can be expressed as linear functions of 3 uncorrelated normally 
distributed random variables, Y x , Y 2 , and Y 3 . And, since uncorrelated normally 


*^By definition cta a ~ P\ \ a \ cr \ where cu a is the covariance between AR 

Ai r o Ax V Q AAf 0 IAV 0 A\r 0 /_W o Alr o AWQ o 

and AV o , is the coefficient of correlation between AR q and AV q , and cr^ f and cr ^ 

are the standard deviations of AR q and AV q respectively. 

(3)\v T 

is the transpose of W, A 1 is the inverse of A and | A | is the determinant of A. 
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distributed random variables are necessarily independent (in the probability 
sense), their joint probability density function is given by the product of their 
marginal distributions. This is the essence of the technique. First, a linear 
transformation is made from the correlated space of AR , Av . and Ar to an 

r o o o 

uncorrelated space of Y 1 , Y 2 , and Y 3 . Computations are made in the uncorre- 
lated space to obtain the joint probability density function of Y x , Y 2 , and Y 3 . 
Then, using the inverse linear transformation, the joint density of AR o , AV o , and 
Ar o is obtained. Finally, the probability density function of perigee error AR 
is calculated using (B.6) and values for r , v , and y . The cumulative distri- 

O O O 

bution function is obtained by summing the probability density function of AR p . 


Let us first express AR , AV , and Ar in terms of standardized normal 

0 °/ A \ O 

random variables X t , X 2 , and X 3 - ( > The reason for doing this is to avoid 
numerical problems in the machine program which diagonalizes the covariance 
matrix. 


= ■'Ar.Xs 


(B.7) 


From (B.7) it is seen that the coefficients of correlation between the vari- 
ables Xj, X 2 , and X 3 are the same as the coefficients of correlation between the 
variables AR o , AV o , and Ar o . For example, the coefficient of correlation between 
X. and X, , p is the same as the coefficient of correlation between AR„ and 

1 / e \ X 1 X 2 ° 

AV , since ( 5 ' 

o * 


p 


Cl X 2 


£ (°a, x. °Av 


X 1 x 2 


a a 

X 1 x 2 


a Av 


T Ar ft ^Av. 


(B.8) 


^^The standardized normal random variable X has zero mean and unit variance, and probability 
density function g(x) given by: 

g(x) = ( 2 tt )~ 1/2 e " x2/2 (-00 < x < 00 ). 

(5) £(Xj) is the expectation of a random variable Xj. £ is the expectation operator. 
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e (ar o av o > 


^At a Av 


(B.8 

cont’d.) 


°Ar„ °Av 


°Ar„Av„ 


°Ar„ °Av„ 


Denoting by X the random vector 


X 


and by A the matrix 


A 



V 


0 



0 


CT Av 0 

0 


0 A 

0 


7 Ay, 


■J 


(B.9) 


(B.10) 


the transformation indicated in (B.7) can be written as 

W = AX (B.ll) 

Thus X has a trivariate normal distribution with zero mean vector and co- 
variance matrix P (also called the correlation matrix) given by (6) 


P = 8(XX T ) = S^' 1 W(A' X W) T ] = £[a _1 WW t (A" 1 ) T ] 

= A" 1 g(WW T )(A- 1 ) T = A” 1 A (A -1 ) T = A -1 AA -1 (B.12) 


^The covariance matrix of a random vector X with zero mean is defined as 8 (XX T ) where 8 


the expectation operator. 
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1 


f 


1 ^Ar 0 Av 0 -°Ar 0 Ay 0 ^ 

^A r „ Av 0 * ^Av 0 A y. 


V 


p kr 0 £v 0 ' 0 Av o Ay o 


A 1 ' 


' 


p A 

X 1 x 2 X x 3 \ 


p 1 p 

^ x l x 2 ^ X 2 X 3 


^ x l x 3 ^x 2 x 3 


v 

(B.12 

cont'd.) 


We now wish to express X as a linear function of a normal random vector 
Y having components which are uncorrelated. Then, since W is a linear function 
of X, the components of W will also be linear functions of these same uncorre- 
lated random variables. Denoting by Y lt Y 2 , and Y 3 uncorrelated normally dis- 
tributed random variables we can write the system of equations as 


*1 

= c n Y i 

+ 

C 12 Y 2 

+ 

C 1 3 Y 3 


X 2 

= <21 Y 1 

+ 

C 2 2 Y 2 

+ 

C 23 Y 3 

(B.13) 

X 3 

= C 31 Y 1 

+ 

C 32 Y 2 

+ 

C 33 Y 3 



or in matrix form 


where 


X = CY 



and C is the nonsingular matrix of coefficients 



r" 

C 1 2 

C 1 3 

C = 

1 

C 2 1 

C 22 

C 23 


\ C 3 1 

C 32 

C 33 


A 

J 


(B.14) 


(B.15) 


(B.16) 
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Since C is a nonsingular square matrix, the inverse B = C" 1 exists, and we 
may write the Y vector as (7) 


Y = C" 1 X = BX 


(B.17) 


Y 1 b ll X 1 + b 12 X 2 + b 13 X 3 

Y 2 = b 21 X l + b 22 X 2 + b 23 X 3 ( B -18) 

Y 3 = b 3 1 X 1 + b 3 2 X 2 + b 3 3 X 3 


Since Y 1 , Y 2 , and Y 3 are uncorrelated, the following relations hold among the 
elements of B and the variances and covariances of X x , X 2 , X 3 and Y x , Y 2 , and 


a 2 
yi 

H b ii 

+ b 12 

+ b l 2 3 + 

2 ( b ll b 12 /0 x 1 x 2 

+ b 1 2 b 1 3 2 x 3 + b 1 1 b 13^ Xl x 3 ) 

a 2 

y2 

= b 2 2 l 

+ h 22 

+ b 2 2 3 + 

2 (b 2 x b 2 2 p xi X2 

+ b 22 b 23 /°x 2 x 3 + b 21 b 23 p> XiX3 ) 

or 2 
V3 

= b 31 

+ b 3 2 2 

+ b 33 + 

2 ( b 31 b 32 ^x lX2 

+ b 32 b 33 ^x 2 x 3 + b 31 b 33^ Xl x 3 ) 


= b n b 21 + b 12 b 22 + b 13 b 23 + (bjj b 22 + b 21 b 12 ) p xiX2 (B. 19) 

+ (t> 12 b 23 + b 22 b 13 ) Px 2 x 3 + ( b ll b 23 + b 2 1 b 13 ) ^Xjxj ® 

b 2 1 b 3 1 + b 2 2 b 32 + b 23 b 33 + ( b 21 b 32 + b 3 1 b 22^ ^x,x 2 

+ (b 22 b 33 + b 32 b 23 ) Px 2 x 3 + ( b 21 b 3 3 + b 3 1 b 23 ) ^ Xl x 3 " ^ 


is the element of B in the i 1 * 1 row and j th 


column. 
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r y,y 3 b H b 31 + b 12 b 32 + b 13 b 33 + ( b U b 32 + b 3 1 b l 2 > Py . , x 2 ( B .19 

cont'd.) 

+ ( b 12 b 33 + b 32 b 13)^x 2 x 3 + ( b U b 33 + b 31 b 13)^ Xl x 3 = 0 

In matrix notation (B.19) is equivalent to the requirement that 

A A 


Y - 


have a covariance matrix Q given by 


Q = 8(YY T ) = 8 


Y Y 
1 2 *1 


V 


Y Y 

3 1 1 


y x y 2 


Y 2 
*2 


y 3 y 2 


r iA\ 


Y Y 

1 2 *3 


Y 2 

I 3 


J 


(B.20) 



^ £(v, 2 ) 

g^Yj) 

8 (Y 

1 Y 3>^ 



= 

e <v,y 2 ) 

8 ( y 2 2 ) 

8 (Y 

2 Y 3 ) 




^ (V, V,) 

s (Y 2 y 3 ) 

e (Y 3 2 ) y 




(< 

a yiy 2 ^yiysA 


(< 

0 

°\ 

— 

a 

yiV2 

cr 2 cr 

y 2 y2y3 

= 

0 

cr 2 

V2 

0 


\ cr 

y yiya 

cr a 2 I 

y 2 y3 y 3 / 



0 

Cr 2 j 

y y 


(B.21) 
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However 0 can be written also as 


Q = 8(YY t ) = B£(XX t )B t = CT'P(CT'f 


(B.22) 


Thus, the problem becomes that of finding a nonsingular matrix which 
diagonalizes P. From matrix theory, since P is a real symmetric square 
matrix of order 3, there exists an orthogonal matrix E of order 3 such that 
EPE t = D, where 




D = 


■ 


V 


is a diagonal matrix of order 3. 




A 




(B.23) 


j 


Moreover, since P is also positive definite, the diagonal elements (eigen- 
values) X- j, X 2 , and X 3 are all positive. Any diagonal matrix which is obtained 
by another orthogonal transformation of P, has the same diagonal elements \ v 
A- 2 , and , possibly in a different arrangement. Hence, in the above we can 
let B be an orthogonal matrix. Then \ x , \ 2 , and \ 3 are the variances of Y x , 
Y 2 , and Y 3 ,ora y 2 i , a y \, a y \. 


As an example, let P AfoAvo - P Avo A ro 
P is given by 

A p 




AAr 0 Ay 0 

A 


p. (0 < p < 1) Then 


(B.24) 


'J 


To find the eigenvalues of P we must find the roots of the secular equation 
(characteristic equation) 


P - Ml = 0 


(B.25) 
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Here I is the identity matrix of order 3. Note that the eigenvalues of a 
real symmetric matrix are always real. 

Equation (B.25) can be written in expanded form as 

(\-l) 3 - 2 p 3 - 3p 2 (X- 1) = 0 (B.26) 


Solving (B.26) for the roots we obtain 

k l = a y] ~ 1 + 2 p 

^ ^ = 1-P 

\ 3 = cr y 2 = 1 - p 


(B.27) 


Since two of the roots are equal, the orthogonal matrix for transforming P 
into diagonal form is not unique. However, it can be constructed. Let us now 
construct the orthogonal matrix C. Since C is orthogonal, C" 1 = C T . Therefore, 
C = B t or 



'21 



12 



c 


'12 


'22 


J 3 2 


21 


'22 


'23 


(B.28) 



'23 


33 , 



'32 



The eigenvector, x x = (c u c 21 c 31 ) T corresponding to = (1 + 2 p) must 
be such that the components c x (j = 1, 2, 3) must satisfy the matrix equation 


PXj “ 
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or 


A p A 

P I/O 

V ' 7 

or expanding the above 

- 2c u + 


A 



1 

C 2 1 

= (l+2/o) 

C 2 1 

W 


W 

C 2 1 + 

C 3 1 = 0 



1 1 


- 2c 


2 1 


'3 1 


' 1 1 


'2 1 


- 2 c 


3 1 


= 0 


(B.29) 


(B.30) 


From the first two equations of (B.30) we obtain c lt = c 21 and from the 
second two equations c 21 = c 31 . Since equation (B.29) is linear, we may make 
the length of the vector Xj unity, thus requiring 

c n + c 2 i + 'si = 1 < B - 3 » 

Equation (B.31) is satisfied if 


c 1 1 c 2 1 c 3 1 



(B.32) 


The eigenvector 
components of c j 2 (j 


x 2 corresponding to \ 2 = (1 - p) must be such that the 
= 1, 2, 3) satisfy the matrix equation 


f c '\ 

l 

( ' z ”\ 

1 

C 2 2 

ii 

/-"N 

l 

C 2 2 

V° 3! 7 

i 1 

W 


(B.33) 
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from which we obtain three equivalent scalar equations in the c ; 2 ’s 


C 1 2 

+ C 2 2 + C 3 2 = 

0 

C 1 2 

+ C 2 2 + C 3 2 = 

0 

C 1 2 

+ C 2 2 + C 3 2 = 

0 


^12^ 


a 13 a 

We may choose for x 2 = 

C 2 2 

and x 3 = 

C 23 


^32 J 


V 33 J 


(B.34) 


any two linearly inde- 


pendent vectors subject to the condition that they are orthogonal and that they 
satisfy the above equations in (B.34). Setting c 12 = 0 in (B.34) we obtain c 22 = 


-c 32 . Thus, if c 22 = y=- and c 32 = - ^=, x 2 will be normalized, 
we will have the dot product (x t ) • (x 2 ) = 0.^ 8) Hence 


Furthermore, 


f°\ 

% = (B.35) 

w 

To obtain the components of x 3 we use the fact that we must have x 1 • x 3 = 0 
and x 2 • x 3 = 0 simultaneously, that is 


C 1 3 + C 2 3 + C 33 


= 0 


C 23 C 3 3 


(B.36) 


should be noted that the eigenvectors of a real symmetric matrix are necessarily orthogonal. 
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from which we obtain 


C 1 3 " _ 2 C 2 3 


and 


C 2 3 C 3 3 


Normalizing we obtain 



Vvey 


Thus, the orthogonal matrix C is 


C 


and therefore B is 


B 


f± 

0 

-A 

vs 


ft 

1 

1 

i 

"ft 

ft 

ft 

1 

1 

i / 


~~ft 

ft) 


Hr 

1 

~ft 

A 

0 

1 

i 

w 

~w 

2 

1 

1 ; 

\> 

w 

wj 


(B.37) 


(B.38) 


(B.39) 


(B.40) 


It can be verified that the elements of B satisfy the equations in (B.19) which 
relate the variances and covariances between the components of X and Y with 
P X[I2 = /°x 2 x 3 - P* lS j = p. The linear equations for this example are there- 
fore, from (B.13) 
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(B.41) 


(- w) ¥ ‘ + ( °> y = + {-wh 

fw) ¥ ' t (w) Y ^(^) ¥ 3 

fwH* (-w) Y * + (ir) Y < 


^x I x 2 * ^2*3 


For the other combinations of signed values for the correlation coefficients 
, and /0 XiX3 which keep the covariance matrix of X nonsingular we 

have the following orthogonal matrices Cj , C 2 , and C 3 which diagonalize P t , 

P 2 , and P 3 respectively. 


For the covariance matrix 

f 1 

■P 

K p 

the orthogonal matrix 


■p 

1 

P 


-A 


!/ 


V5 

0 

' ye 

1 

1 

1 

ya 

" i2 

■ ye 

1 

1 

1 

w 

/ 2 

" ye . 


(B.42) 


(B.43) 


diagonalizes P x . 

For the covariance matrix 

f 


1 
~ P 

\ P 


■p 

1 


A 


V 


(B.44) 
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the orthogonal matrix 



f 

1 

V3 

o 

a 

' i6 

\ 



1 

1 

1 


(B.45) 

C 2 = 

i 

V3 

72 

~7T 




1 

1 

i 



V 

77 

Y2 

Wj 

i 

diagonalizes P 2 . 






And, for the covariance matrix 






f 1 

P 

-A 


P 3 = 


p 

1 

- p 

(B.46) 



\ : p 

~ P 

V 


the orthogonal matrix 







f- 

i 

w 

o 

■A 




i 

1 

1 

(B.47) 

c 3 = 

"■ 

V3" 

~7T 

ft 


l , 

1 

i 

1 J 




w 

"7 2 




diagonalizes P 3 . 

Each of the orthogonal matrices indicated above was obtained in the same 
manner. It should be noted that in this special case where each coefficient of 
correlation has the same absolute value, each covariance matrix has the same 
eigenvalues: \ 1 = 1 + 2p, k 2 = 1 - p , \ 3 = 1 - p, inasmuch as the char- 

acteristic equation is the same in each instance. 

In the same manner as was done for the random vector W we may express 
Y in terms of a standardized normal vector Z having uncorrelated standardized 
normal random variables Zj, Z 2 , and Z 3 as components. 


54 



(B.48) 


Y, = cr Z, 
1 y 1 1 


Y 2 °y 2 Z 2 


Y 3 CT y 3 Z 3 


Finally, we may write AR o , AV o , and AF o as linear combinations of Z 2 , Z 2 , 
and Z 3 . 


A R 0 

AV 

O 

Ar 

o 


( CT Ar 0 CT yi ) C U Z l + 

( a Av 0 CT yi ) C 2 1 Z l + 

( a A ro CT yi ) c 3i Z l + 


( a Ar 0 CT y2 ) C 12 Z 2 
( CT Av 0 CT y2 ) C 2 2 Z 2 
( a &y 0 a y 2 ) c 3 2 Z 2 


+ Kr 0 CT y3 ) C 13 Z 3 
+ K/ys) C 23 Z 3 
+ (° r A ro CT y3 ) c 3 3 Z 3 


(B.49) 


Now, since Y, , Y„ , and Y, are uncorrelated normal random variables, they 

12 o q r _ 

are independent in the probability sense. ( ' Therefore, their joint probability 
density function f (y 3 , y 2 , y 3 ) factors into the product of the marginal distribu- 
tions of Y X ,Y 2 , and Y 3 , or f ^y^, f 2 (y 2 ), and f 3 (y 3 ). 

f(yi> y 2 *y 3 ) = f i(yi) f 2 (Y2) f 3(y 3 ) (B.50) 

Denoting by g (x 3 , x 2 , x 3 ) the density of X 3 , X 2 , X 3 , and noting that by (B.13) 
the transformation from the y -space to the x-space is one-to-one, as well as 
from the x-space to the y -space (B.18), the density of X,, X 2 and X 3 is 

g(x 1 ,x 2 ,x 3 ) = f[y 1 (x 1 ,x 2 ,x 3 ),y 2 (x 1 ,x 2 ,x 3 ),y 3 (x 1 ,x 2 ,x 3 )]j(x 1 x 2 ,x 3 ) (B.51) 


(9) 


We note that uncorrelated random variables are not always independent. 


random variables are always uncorrelated. 


However, independent 
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where J (Xj.X 2 .X 3 ) is the Jacobian 


J (xj, x 2 , x 3 ) - mod 


111 

B Yi 

hi 

3 x j 

d x 2 

dx 3 

111 

^y 2 

hi 

d x j 

d x 2 

B x 3 

b y 3 

hi 


3 x j 

B x 2 

B x 3 


(B.52) 


and "mod" means the absolute value of the expression following it. Using (B.18) 
we can compute J (Xj, x 2 , x 3 ) as follows: 


J ( x l« x 2 > X 3 > 


mod 


1 1 

b l 2 

b 13 



2 1 

b 22 

b 23 

mod | B | 

(B.53) 

3 1 

b 32 

b 33 




However, since B is an orthogonal matrix, J (Xj, x 2 , x 3 ) = 1. Therefore, 
using (B.50) and (B.51) the joint density of X 1# X 2 , and X 3 with the aid of (B.13) 
is given by 

P{xj < Xj £ Xj + Axj, X 2 < x 2 < X 2 + Ax 2 , x 3 < X 3 £ x 3 + Ax 3 } 

= g(Xj,x 2 ,x 3 )AxjAx 2 Ax 3 = f (yj,y 2 ,y 3 ) AyjAy 2 Ay 3 

= P{y t < Yj < y t + Ay j , y 2 < Y 2 < y 2 + Ay 2 ,y 3 < Y 3 < y 3 + Ay 3 } Ay j Ay 2 Ay 3 

= P{cj jXi + c 2 jX 2 + c 3 iX 3 <Yi<Cji(Xj + Axj) + c 2 j(x 2 +Ax 2 ) + c 31 (x 3 +Ax 3 )} 

* P { C 12 X l + C 22 X 2 + C 32 X 3 <Y 2- C 12( X l +Ax l) + C 22( X 2 +Ax 2) + C 32( X 3 +Ax 3>} 

X P{ C 13 X 1 + C 23 X 2 + C 33 X 3 <Y 3 — C 13( X 1 + Ax l) + C 23( X 2 + Ax 2^ + C 3 3^ X 3 + Ax 3 )} 

x A yi A y 2 A Y3 
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X 


f l( c ll x 


+ C 21 X 2 


+ C 31 X 3)^2^ C 12 X 1 +C 22 X 2 + C 3 2 X 3 
x ( c i l Ax i + c 2 l AjC 2 + C 3 l Ax 3 ) * 


) f 3 ( C 13 X l + c 23 X 2 + C 33 X 3) 
( c 12 Axj + c 22 Ax 2 + c 3 2 Ax 3 ) 


( C 13 AX 1 + C 23 AX 2 + C 33 AX 3 ) 


(B.54) 


For purposes of machine computation, the range of each random variable 
Yj (i = 1, 2, 3) was approximated by (2N + 1) mutually exclusive intervals of 

length (±)a yi extending from -Lo-^ to Lo- yj , and the probabilities of falling 

within these intervals approximated in the following manner. 


First, (2N + 1) values for the normal probability density function were cal- 
culated for each Y { 


where 


= (t) t k -< N * »] 

(k = 1,2,3, .... ( 2N + 1)) 
(i = 1,2,3) 


Then, the sum S. was formed 

7 i 


2N+ 1 



k = 1 


(B.55) 


(B.56) 


Next, the q. 's were normalized, i.e., p 's were formed where 

1 k 


lk a i (B.57) 

(k = 1,2,3 (2N + 1)) 

(i = 1,2,3) 
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Next, defining Y' x , Y' 2 , and Yj as discrete random variables approximating 
the normal random variables Yj , Y 2 , and Y 3 we may write 

p ( y : = *;„) = < B - 58) 

(i - 1,2,3) 

(k = 1,2,3 ( 2N + 1)) 

And, we may write the probability that Y i should fall within the interval 

(*:„ ♦(yK) (B-59) 


where 


Ay, 



as follows 


p f: t -(4K, + (tK}- p ( y ; = 

(i = 1,2,3) 

(k = 1, 2, 3, . . . , ( 2N + 1)) 


(B.60) 


Thus, each normally distributed random variable Y i (i = 1,2,3) has been 
approximated by a discrete random variable Y? (i = 1, 2, 3) in the range -L a 

to +Lcr and having probability mass points p. (i = l,2,3;k=l,2,3,... 

y i lk l /1\ 

(2 N + 1)) at the center of each subdivision (yj - ^-y^Ay . , + (yjAy.). 

It is easily verified that the probabilities for each sum to 1, since 


2N+ 1 


2N+ 1 


2N+ 1 


L p{y: = = Z>„ = (<r)ZS 

k= 1 k= 1 k= 1 


(B.61) 
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Since each Y! can assume (2N + 1) possible values, y| t , each with probability 

P ifc (k = 1, 2, . . . 2N + 1) there will be a total of (2N + l) 3 possible values for 

the approximation to the joint probability density function of Y t , Y 2 , and Y 3 . 

Since the Y. 's are independent, the probability of occurrence k { of each com- 
bination of values for yj , yj , and yj (j ,k, 1 = 1, 2, 3, . . ., (2N + 1)) is 

given by the product of p, , p~ , and p, , i.e., 

A j z k 

Pjkl “ Plj P 2 k P3, (B.62) 

(j.k.l = 1,2,3, .... ( 2N + 1)) 

In order to distinguish between the continuous random variables and the dis- 
crete approximations to these variables, let Xj, Xj, Xj, AR^, A\P, AP, and AR p 
be the discrete random variables which approximate x 1 ,x 2 ,x 3 ,ar o , av o , Ar o , 
and AR p , respectively. 

For each combination of values for Y| , YJ, and Yj , the probability of occur- 
rence was obtained using (B.62), the inverse transformation to the correlated 
set XJ , X' 2 , and Xj made using (B.13), and values calculated for AR^, A\P, and 
AT p using (B.7). Since there is a one-to-one correspondence from the uncorre- 
lated y -space to the correlated space of AR^, A\P, and AP, the joint density 
function approximation of AR q , AV o , and Af o will have (2N + l) 3 possible values, 
and for each combination of values for AR^, A\P , and AP^ we may obtain a value 
for AR p , using (B.6). Therefore, the discrete approximation to the probability 
density function of AR p will also have (2N + l) 3 possible values. At this point 
we should note that AR p as given in (B.6) depends on actual values r o , v o , and 
y • 


We will now discuss how the actual values for r o , v o , and y o were calcu- 
lated. 

First, actual values for the semi -major axis a, the true anomaly of the 
spacecraft at insertion B , and the eccentricity of the parking orbit e , were 
assumed. 

The insertion radius r o was obtained using (A.3) of Appendix A. 


r 


O 


a ( 1 - e 2 ) 

1 + e cos 8 


(B.63) 
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The insertion speed v o was obtained from (A.l) of Appendix A. 



(B.64) 


In order to obtain the flight path angle at insertion y o , let us rewrite (A. 2) 
of Appendix A. 


1 - 



From (A.l) 


r 

O 





(B.65) 


(B.66) 


Substituting (B.66) into (B.65) gives 


1 - e 2 = 


-3G) 


cos y n 


(B.67) 


Solving for 



in (B.63) and substituting the result in (B.67), we have 


/ e sin 6 \ 

= tan ri — cosTi 


(B.68) 


The actual (or assumed) perigee radius r p was obtained using (B.l). 

Thus, with assumed actual values for r Q > v 0 » 7 0 > and r p , and by using all 
the (2N + l) 3 possible combinations of values for Yj , Y 2 , and Yj with their 
corresponding (2N + l) 3 probabilities p jkl = p 2k P 3l (j » k , 1 = 1, 2, 3, . . ., 

(2N + 1)), an approximation to the probability density function of AR p through 
the transformations indicated in (B.13), (B.7), and (B.6) was obtained. 
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It is easily seen that the probability density function approximation for AR p 
sums to 1 . 


2 IN + 1 

f (Ar p ) d(Ar p ) - 2^ Pijk 

i » j t k — 1 


2N+ 1 



i , j , k- 1 







(B.69) 


where f (Ar p ) is the probability density function of AR p . 

The cumulative distribution function of AR p orF (Ar p ) is given by 

Ar p 

f (x) dx (B.70) 


F(Ar ) = 




For purposes of machine computation, the range of AR p was subdivided into 
2000 mutually exclusive intervals or "bins" one tenth of a kilometer wide and 
extending from -80 kilometers to +120 kilometers. For all values of AR^ falling 
within one of these intervals, the associated probabilities were summed. This 
can be done since all of the (2N + l) 3 possible combinations of values for AR p , 
AV p , and AT p are mutually exclusive. In this fashion we obtain an approximation 
to the probability density function of AR p . The cumulative distribution function 
approximation for AR p was then obtained, by summing the probability density 
function approximation for AR p . 

A Fortran IV computer program was written and used to calculate the 
probability distribution of perigee error on the UNIVAC 1108 computer. For 
L = 4 and N = 50, i.e., approximating the normal curve from -4 standard de- 
viations to +4 standard deviations in (2N + 1) = 101 equally spaced intervals 
(0.08 standard deviations in width), the calculations took approximately 4 minutes 
to obtain the probability density function and cumulative distribution function ap- 
proximation to AR p for a given set of insertion conditions. 
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Exp re ssion For the Error in Perigee For a Circular Orbit 


The expression for the error in perigee for a circular orbit is somewhat 
less complicated than the expression in (B.2). We will now treat the circular 
case. 

Let us first rewrite equations (A.l), (A.2), and (A.4) of Appendix A. 



r - a ( 1 ~ e) . 

p v ' 

Substitute (B.72) into (B.71) and also into (B.73) 


1 - e 2 


v 2 r 2 cos 2 y 

o o 7 o 



(B.71) 


(B.72) 


(B.73) 


(B.74) 


r 

p 



1 - e 


(B.75) 


or 


e - 



(B.76) 


Rewrite (B.34) and multiply by r . 

( ” u 2 r „ 2 cos 2 7 „ \ j 2 V \ 

r„(l - e)(l + e) = [ ~ j S ^ (B.77) 
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Cancel (1 - e) in (B.77) using (B.75), and substitute (B.76) into (B.77) for 
(e + 1). The result is 


2r , + ,AAAAL) = 

° P \ jJL ) p 




(B.78) 


Substituting ( r o + Ar o ) for r o , ( Vq + Av o ) for Vq , ( y o + Ay^ for y o . 


(r p + Ar p ) for r p , and collecting terms, we have 


2r + [ — ) v 2 r r - 4 r Ar + 

1° \ fX j ° ° P P P 


}ik 


2 + 


v r 

o p 




jJL J ° P 


A v A r + 

0 p 


Ar* Ar 

° p 


(-M. * r - 2 ' 

\ jJL / O O 


(Ar p ) : 


” 2 r - f ) v 2 r 2 + ( ) v 2 r 2 cos 2 

P \/X/° P \ fl / o o 


y 

7 n 


Ar 


(ih 


, . r 3 cos 2 y 

/Li J 00 7 o 


-(2-)v r r 2 

\ jj. / O o p 


Av 


[■ (t) v . 


2 r 3 sin y cos 7 Ay 

> o / o 7 o 7 0 


o 




V r 1 J r 

o p 


(JJ\ 2 2 

) v r z cos z y 

\ UL I O O 7 0 


Ar Av 

o o 


[' (7) v 


r 3 sin y n cos y Av Ay 

00 7 o 7 o o 7 o 


3 


+ 



cos y 

7 o 


A r A y 

o 7 o 
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The minus sign has to be used in (B.81) since Ay o has to produce a negative 
Ar p if Ar o = Av o = 0. (The plus sign is applicable to the error in apogee.) 
Thus, the perigee error for a circular orbit becomes: 
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Derivation of the Probability Density Function of Perigee Error For a 
Circular Parking Orbit When the Insertion Errors are Perfectly Corre- 
lated and Normally Distributed. 


As a check against the calculations performed by the computer, the prob- 
ability density function of the error in perigee for a circular orbit was obtained, 
assuming the insertion errors to be perfectly correlated (coefficients of corre- 
lation between the insertion errors of +1.0) and normally distributed. 


If AR q , AV o and AF o are normally distributed random variables mutually 
correlated with coefficients of correlation of +1.0 , there is complete linear de- 
pendence between them. Thus, they will vary in the same sense. Let X denote 
the standardized normal random variable (zero mean and unit variance). We 
may then write (B.82), substituting AR p for Ar p , AR o for Ar o ,AV o for Av o , and 
♦Ar o for A y o , as 


AR p 


2 



+ AR 


•] 



(r Ar ) 

v o o ' 


2 



(B.83) 

where a A , o■^ , cr* are the standard deviations of AR , AV and AF respec- 

Ar 0 Av 0 ’ /xy a o o o 

tively. 


Letting 1 



and 


m 



< r o°A? 0 ) 
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we can write (B.83) as 


AR p = IX - m | X| (B.84) 

Letting F(Ar p ) denote the distribution function of AR P , G(x) the distribution 
function of X, and assuming m > 1, we have: 

F(Ar p ) = Prob{AR p < Ar p } = Prob{l X - m |X | < Ar p } (B.85) 

For X < 0, |x| = -X, and we have 

Prob{lX-m|x| < Ar p } = Probjx < yyj = G (jy) (B.86) 

For X 2 0, |x| = X, and we have 
Prob{l X - m | X| < Ar p } = Prob{(l -m)X<Ar p } 


- Prob •[(m - l)X>-Ar p } 


1 - Prob 



- Ar 

P_ 

(m - 1) 


} 


= 1 - G 


- A: 


(m - 1) 


(B.87) 


Since the events X < o and X > o are mutually exclusive, the total probability, 
or F (A r p ) is given by the sum of (B.86) and (B.87) 


F(Ar ) = Prob{lX - m |X| < Ar } 


- Prob 


>{ Xi TT^)} + Prob { x j T^T)} 


(B.88) 
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PROBABILITY / (km) 


(B.88) 

(cont'd.) 


= 1 + G 


t *r_ \ / -ir. \ 


Letting f (Ar p ) be the probability density function of AR p> and g(x) the 
standardized normal density function, then we have the following relationship 
between f(Ar p )and g(x). 


f(Ar p ) = F'(Ar p ) 



(B.89) 


These functions f (A r p ) and F(Ar p ) are shown in Figures 17 and 18. 


f (Ar p ) 



(km) 

Figure 17. The Probability Density Function f of Perigee Error AR p for a Circular Orbit when 
the Insertion Errors Are Perfectly Correlated ( /0 Ar o Av o = PAv Ay = PA t A y ~ +1.0). 
Insertion Parameters: a = 3544 n.mi. (6563 km), e = 0. 
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1.0 



Figure 18. The Cumulative Distribution Function F of Perigee Error AR p for a Circular Orbit 
when the Insertion Errors Are Perfectly Correlated (£>A r Av “ ^°Av Ay = /°Ar Ay “ 1*0)* 

o o o o o ' o 

Insertion Parameters: a - 3544 n.mi. (6563 km), e = 0.0. 


Expression for the Perigee Error for Eccentricities In the Range 
0.005 S. e < < 0.05 


For eccentricities in the range of approximately 0.005 to 0.05, the following 
expression taken from Reference 5 (equation (9)) is a good approximation for the 
error in perigee. It is repeated here. 


AR % 

p 


(2 - cos $)AR q + 2(1 - cos 



- (r o sin #)Ar o 


(B.90) 


It is evident that, if AR 0 , AV oJ and Ar o are all normally distributed, then, for 
a fixed 6 (true anomaly of the spacecraft at insertion) and e (eccentricity), 
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AR p will also be normally distributed, having a zero mean, and a variance given 

by 


o\^ r - (2 - cos 9 ) 2 <j£ r + 4(1 - cos 9) 


/ r V 
I ° 1 

\ v / 

\ ° / 


°a * 0 


+ r Q 2 sin 2 9 cr£, + 4(2 - cos 9) (1 - cos 9) 




0 Av 0 ^Aro a Av 0 


-4(1- cos 9 ) 



sin0) ^Av 0 A yo °-Av 0 ^A ro 


(B.91) 


- 2(2 - cos 9) (r o sin 5) PAr 0 Ar 0 °Ar 0 CT A 7o 


In order to compare the distribution of AR p as calculated by numerical in- 
tegration with the normal distribution approximation, we note Figures 19 and 20. 
In each instance, the semi-major axis is a = 6563 km, the eccentricity e =0.01, 
and the standard deviations of AR o , AV 0 , and Ar o are 

cr^ r = 0.8 n.mi. (1.4816 km) 


CT A V 


O 



ft /second (1.62 m/sec) 


(B.92) 


a A 


To 



(9.3 mrad) 


In Figure 19, the coefficients of correlation between the insertion errors 
are all +0.9. In Figure 20, the insertion errors are uncorrelated. Both figures 
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show the 99.5% probability points as a function of true anomaly 8 . The agree- 
ment between the two curves is good. 



Figure 19. Comparison Between the Required 99.5% Padding as Calculated by Numerical 
Integration and the Required Padding Using a Normal Distribution Approximation. Insertion 
Parameters: a = 3544 n. mi. (6563 km), e = 0.01. Insertion Errors Are Positively Correlated 

(Coefficients of Correlation = +0.9). 
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Figure 20. Comparison Between the Required 99.5% Padding as Calculated by Numerical 
Integration and the Required Padding Using a Normal Distribution Approximation. Insertion 
Parameters: a - 3544 n.mi. (6563 km), e - 0.01. Insertion Errors Are Uncorrelated. 
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APPENDIX C 


THE CUMULATIVE DISTRIBUTION FUNCTION OF THE ERROR IN 
ECCENTRICITY. 

The eccentricity e can be expressed as a function of r o , v o , and y o . From 
Appendix A (equation A. 6) 


e - e(r ol v q> y c ) 



(C.l) 


In calculating the insertion parameters r o , v o , y 0 , the errors Ar o , Av o , and 
Ay o are introduced. Therefore, the calculated eccentricity e ca x will deviate 
from the actual eccentricity e , by the amount Ae 

Ae = e cal ~ e = e ( r o + Ar o’ V o + Av o’ + ^o) " e ( r o- V o> T 0 ) (C.2) 


Letting AE be a random variable representing the error in eccentricity, and 
AR o , AV o , AT o random variables representing the errors in insertion radius, 
speed, and flight path angle respectively, we may write with the aid of (C.2) 


AE - e (r + AR 

v O O 


+ AV , y + AF ) “ e (r 

O ’ ' O O 7 V f 


T 0 ) 


(C.3) 


The probability density function approximation for AE was obtained in an 
analogous manner as was done for AR (Appendix B). In an effort to shorten 
machine calculation time, however, each normal probability density function 
(in the uncorrelated y -space) was approximated by a discrete probability dis- 
tribution having 81 mass points (instead of 101), extending from -4 standard 
deviations to +4 standard deviations. The width of each bin for the error in ec- 
centricity was taken as 10 " 5 . 

The Error In Eccentricity for a Circular Orbit 


Equation (C.l) may be written 


1 - e 2 


2r v 2 

O O 


cos" y 



(C.4) 
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Substitute (e + Ae) for e . ( r + Ar ) for r (\r + Av i for 

' ' ' ' O O' ~ O * v ’ O of 

for y 0 and collect terms. 


\T 

O ’ 


/ o/ 

\'o 


/W 


\ 


1 - e 2 - 2eAe - (Ae) 2 



cos 2 y 

' o 


+ 





cos (2 y o ) 


(^y 0 ) 2 



sin y cos y 

' O ' o 



sin y 

' o 


cos y 


A r c 
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r v 


o o 


cos 2 y 

' o 
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cos 



Av 

o 



sin y cos y 

' O ' o 



sin y cos y 

7 O / o 


Av Ay 

o / o 





Av 2 


O 


+ {(“yH 2cos27 ° r o v o 4 c ° s 2 ^oj Ar o 


+ 



r v 4 sin y cosy 

o o / o ' o 



sin y Q cos y 


Ar Ay 

o ' o 



CO S 



Av 
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-ft) 


2 4 2 

r v co s y 

° ° ° (C.5) 
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1 - (Ae ) 2 = 2 - (A 7q ) 2 (Av o ) 2 - Ar o Av o - 1 (C.6) 



O 


Solving for Ae in (C.6), we have, 



Since e = 0 for a circular orbit and the calculated eccentricity must be 
positive (by the definition of eccentricity), Ae in (C.7) must also be positive. 
Therefore, the error in eccentricity is, 
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APPENDIX D 


THE CUMULATIVE DISTRIBUTION FUNCTION OF THE ERROR IN 
TRUE ANOMALY. 

Using (A.3) we may calculate 6 


e - 



(D.l) 


where a and e are functions of r o , v o , andy o by (A. 5) and (A. 6) 



(A.5) 


e 



(A.6) 


The correct quadrant for 6 may be resolved from inspection of (B.68) 


tan y 

' o 


e sin 6 
1 + e cos 6 


(B.68) 


We see that for y o > 0°, 0 £ 6 < 180° and for y Q < 0°, -180° < 0 < 0°. 
Thus, we may write 


d ~ d (r D , v o , y 0 ) (D.2) 

Letting AR o , AV o , IT be random variables representing the errors in in- 
sertion radius, speed, and flight path angle, with the aid of (D.2) we may write 

A0 - 6 (r +AR,v + A V , y + Ar ) - 6 (r , v , y ) (D.3) 

v o O O O 'O O' v O O 'O' 
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where A0 is a random variable representing the error in true anomaly for a 
non-circular orbit. 

The cumulative distribution function of A0 

K (A0) = P (A0 < A0) (D.4) 

may be seen in Figures 21 and 22. These were obtained in analogous fashion as 
for AR p and AE . 

Figure 21 is for uncorrelated insertion errors, an eccentricity of e = 0.001 
and true anomalies at insertion of 6 = 0°, 90°, 180°, and 270°. 

F igures 22a and 22b are for the same conditions of true anomaly at inser- 
tion and correlation, except the eccentricity is 0.005. 



Figure 21. Cumulative Distribution Functions of the Error in True Anomaly for Various Actual 
Values of True Anomaly at Insertion. Insertion Parameters: a = 3544 n.mi. (6563 km), e = 

0.001. Uncorrelated Insertion Errors. 
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(degrees) 

Figure 22a. Cumulative Distribution Functions of the Error in True Anomaly for Various Actual 
Values of True Anomaly at Insertion. Insertion Parameters: a = 3544 n.mi. (6563 km), e ” 

0.005. Uncorrelated Insertion Errors. 



Ar 0 Ay 0 =0.9 


A 0 (degrees) 

Figure 22b. Cumulative Distribution Functions of the Error in True Anomaly for Various Actual 
Values of True Anomaly at Insertion. Insertion Parameters: a - 3544 n.mi. (6563 km), e - 

0.005. Correlated Insertion Errors (p& r = +0. 9 ' ^Av Ay ~ P& A y “ ±0 - 9 ) 
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The "3-sigma" limits used for AR q , AV o , and Af o are 



= 2.4 n. mi. (4.44 km) 


3 a Av 0 

= 16 ft/sec (4.87 m/sec) 

(D.5) 

3 °Ar 0 

- 0. 16° (2.79 mrad) 



From inspection of the curves it can be seen that for parking orbits with 
low eccentricities (e < 0.001) true anomaly at insertion 6 cannot be determined 
accurately. For example, due to tracking errors, for e = 0.001, 6 = 0° or 
6 = 180°, and uncorrelated insertion errors (dashed curve of Figure 21) there 
is a 90% probability that the error in true anomaly will be between -66° and 
+ 66 °. 


From Figures 22a and 22b it can be seen that true anomaly can be deter- 
mined more accurately if the eccentricity is greater. 
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